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We have written a Mathematica program that calculates the integrand corresponding to any 
amplitude in the closed-time-path formulation of real time statistical field theory. The program 
is designed so that it can be used by someone with no previous experience with Mathematica. It 
performs the contractions over the tensor indices that appear in real time statistical field theory 
and gives the result in the 1-2, Keldysh or RA basis. We have used the program to calculate the 
\^ . ward identity for the QED 3-point function, the QED 4-point function for two photons and two 

' fermions, and the QED 5-point function for three photons and two fermions. In real time statistical 

\ field theory, there are seven 3-point functions, 15 4-point functions and 31 5-point functions. We 

. produce a table that gives the results for all of these functions. In addition, we give a simple general 

^j^' expression for the KMS conditions between n-point green functions and vertex functions, in both 

r-{ ' the Keldysh and RA bases. 

. PACS numbers: 11.10.Wx,11.15.-q 
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^ ■ I. INTRODUCTION 

OO 

I It is well known that calculations in real time statistical field theory are complicated by the proliferation of indices 
I that results from the doubling of degrees of freedom (for a recent review see [l|). This difficulty causes many people 
' to avoid using real time finite temperature field theory, in spite of its significant advantages over the imaginary time 
\^ , formalism. Two of the major advantages of working in real time are that analytic continuations are not necessary, 
' and that it is easy to generalize to non-equilibrium situations. In this paper we make a contribution towards reducing 
the technical difficulties associated with the real time formulation of statistical field theory. 

We use the closed time path (CTP) formalism of real time statistical field theory which consists of a contour 



Oh! 



I ' with two branches: one runs from minus infinity to infinity along the real axis, the other runs back from infinity to 



Oh, 



minus infinity just below the real axis. This contour results in a doubling of degrees of freedom. Physically, these 
extra contributions come from the additional processes that are present when the system interacts with a medium, 
instead of sitting in a vacuum. As a result of these extra degrees of freedom, n-point functions have a tensor stucture 
which results in calculational complexities that increase geometrically when one considers either calculations at higher 
, loop order, or the calculation of higher n-point functions. Statistical field theory can be formulated in different bases, 
5_] ■ which produce different representations of these tensors. There are three popular bases: the 1-2 basis, the Keldysh 
\ basis, and the RA basis. Most of the original work in this field was done in the 1-2 basis. The Keldysh basis has the 
advantage of being more easily adaptible to non-equilibrium situations. The RA basis produces particularly simple 
expressions in equilibrium. 

This paper is organized as follows: In section^]we review the 1-2 formalism. In section lllll we discuss the issue of ba- 
sis transformations in general. Our discussion follows that of 0|. In sections IVland lvTl we give expressions for n-point 
functions and vertex functions in the Keldysh and RA bases. In section lVIll we give the KMS relations in each of these 
bases. These equations give a set of relations between the various components of a given n-point function or vertex 
function which hold in equilibrium, and are often useful for simplifying the expressions that result after contracting over 
indices. We note that some of these expressions have appeared previously in the literature, using a slightly different 
notation 5,6]. We give general expressions in the Keldysh and RA bases. We emphasize the simplicity of the expres- 
sions in the RA basis. In section rVIIII we discuss the Mathematica computer program that we use to perform most of 
the calculations in this paper. This program is designed so that it can be used by someone with no previous experience 
with Mathematica. It is available on the internet at www.brandonu.ca/physics/fugleberg/Research/Dick.html. 

The program performs contractions over CTP indices and produces the integrand that corresponds to any n-point 
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vertex diagram, in the 1-2, Keldysh or RA basis. Several options to simplify the results are built into the program 
and can be selected by the user as part of the input. For example, for calculations in equilibrium, the KMS conditions 
can be automatically implemented. We discuss the basic design of the program and explain in detail, using a simple 
example, how the program can be used. The program treats all fields as scalar bosons and from the beginning of 
the paper to the end of section IVlIII we discuss only scalar bosons. It is straightforward to apply the results of the 
program to other field theories. In section ITXI we illustrate the use of the program with the calculation of some QED 
ward identities. The general form is the same as the well known zero temperature expression. In real time statistical 
field theory however, the situation is complicated by the additional degrees of freedom. Each n-point function has 
2" — 1 independent components and each of these components has its own ward identity. We start by looking at 
bare 1-loop diagrams and obtaining the form of the ward identities for the 3-point function (2" — 1 = 7 components), 
the 4-point function for two fermions and two photons ((2" — 1 = 15 components), and the 5-point function for two 
fermions and three photons (2" — 1 = 31 components). For the 3- and 4-point functions, we verify that the same 
expressions hold when full corrected vertices and propagators are used. Results for some components of the 3-point 
function were obtained in and . A partial set of expressions for bare 1-loop 3- and 4-point diagrams in the htl 
approximation was derived in "9*1 . In this paper we produce a complete set of expressions. Our results agree with the 
partial results of 0,0,^]; which provides a check of our program. We note that the full calculation involving corrected 
vertices would be prohibitatively difficult without the use of a program like the one presented here. Our conclusions 
are presented in section Ixl 



II. THE 1-2 BASIS 



The ri-point function is defined: 



(1) 



where the subscripts {bi} take values 1 or 2 to indicate which branch of the contour the corresponding time argument 
falls on, and the symbol V represents ordering along the closed time path. In what follows we will suppress the 
superscript (n) and let the number of indices indicate the number of fields in the n-point function. All 7i-point 
functions will be denoted G, except for the 2-point function, or the propagator, which will be called D. 
In the 1-2 basis the 2-point function can be written as a 2 x 2 matrix: 



D(l_2) - 



D21 D22 



(2) 



The component Dn indicates a propagator for fields moving along the top branch of the contour, D12 is the propagator 
for fields moving from the top branch to the bottom branch, etc. In co-ordinate space these components are given by, 

D^^{x,v) = ~i{T{^{x)<P(y))) (3) 

Di2{x,y) = -~i{(j){y)(j){x)) 

D2i{x,y) = ~i{(j>{x)(j>{y)) 

D22{x,y) = -i{f[c^{x)c^{v))) , 

where T is the usual time ordering operator and T is the antichronological time ordering operator. These four 
components satisfy, 



I] = ^11 - D12 - D21 +D22=0 



(4) 



a=l 6=1 



as a consequence of the identity 0{x) + 0{~'x) — 1, and thus only three components are independent. 
The 3-point fimction in the 1-2 basis is a (2 x 2 x 2) tensor with components in co-ordinate space: 



Giii{x,y,z 
Gii2{x,y,z 
Gi2i{x,y,z 
G2ii{x,y,z 
Gi22{x,y,z 



= -(r(0(a;)0(y)0(z))) 

= -{4>{z)T{cf>{x)cf>{y))) 

- ~{cl>{y)T{4>{x)(k{z))) 

= ~{^ix)T{<j>{ymz))) 

= -(r(0(y)0(z))(/.(a;)) 
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G22iix, y,z) = -(r(0(a;)</)(y))0(z)) 

G222ix,y,z) = -{f{(bix)(biy)Hz))). (5) 
Only seven of these components are independent because of the identity 

EEE(-ir'+^+'G.,. = o (6) 

a=l b=l c=l 

which follows in the same way as Q from 9{x) + 6{—x) = 1. 

The 4-point fmiction in the 1-2 basis is a (2 x 2 x 2 x 2) tensor. We write out a few examples: 



Giiii{x, y,z,w 
Giii2{x, y,z,w 
Gii2i{x, y,z,w 
Gi2ii{x,y,z,w 
G2iii{x, y,z,w 
Gii22{x, y,z,w 



t{T<j>{x)<l>{y)^{z)^{w)) 

i{^{w)T{cj,{xmymz))) 
i{cp{z)T{4>{x)4>{y)<k{w))) 
z(0(y)T(0(x)0(0)0M)) 
i(0(a;)r(0(y)0(z)0(«;))) 
z(r(<^(z)<^H)r(0(x)0(2/))) 



(7) 



These functions obey a relation similar to Q and ©: 

Y^Y^Y^Y^i-ir^'^'^'G.^c^-^- (8) 



a=l b=l c=l d=l 

The structure of higher n-point functions is similar. 

Truncated green functions are called vertices and will be denoted T, except for the two point function which will 
be called 11. They are defined by the equation: 

Gbi---b„ — Gij^b^ ■ ■ ■ Gfj^b^r''^ . (9) 



The vertex functions satisfy the constraint: 
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E E-'-Er"^'^ ■'"=0 (10) 



bi = lb2 = l fc„ = l 



III. BASIS TRANSFORMATIONS 



A scattering amplitude is calculated by multiplying together a series of vertices and propagators. By convention 
we will assign lower indices to n-point functions (like the propagator) and upper indices to vertices. Summations are 
carried out over pairs of repeated indices, where one of the indices is an upper index and the other is a lower index. 

The expression for a given scattering amplitude can be transformed to a different basis by performing a rotation. 
There are two commonly used bases: the Keldysh basis and the R/A basis. Both of these bases express results 
in terms of combinations of the components of the 1-2 propagator that have a direct physical interpretation: the 
retarded, advanced and symmetric propagators. These expressions are easy to obtain in co-ordinate space. Using (O 
and I0J one can show: 

Dret{x,y) = Dii{x,y) - Di2{x,y) = -i0{xo - yo){[<l>{x), (f){y)]) (11) 
Dadv{x,y) Dn{x,y) - D2i{x,y) = -i9{yo - Xo){[(l){x),4i{y)]) 
Dsym{x,y) = 1311(2;, y) + D22{x,y) = -«({0(a;), 0(y)) . 

One advantage of the Keldysh basis is that it is easily generalizable to non-equilibrium situations. In equilibrium, 
the R/A basis produces particularly simple expressions for amplitudes, and for expressions like the KMS conditions, 
which give relationships between different amplitudes. 
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The rotation to a different basis is accomplished by matrix multiphcation. We rotate lower indices by multiplying 
by a matrix U and upper indices by multiplying by the matrix V. These matrices are related through the following 
equation. 

Vik) = iU^)-\-k). (12) 

We obtain expressions of the form: 

U-,'^ik)U-,\-k)Dai,{k) - D'^ -,ik) (13) 

Vlikl)V\{k2)Vl{-k, - fc2)r'''"(fcl, fc2, -fcl - fc2) ^ {T'f^'''{ki,k2, -fci - k2) 



where we have simplified the notation by writing a two point function of the form Dab{k, ~k) as Dab{k). It is 
straightforward to see that amplitudes have the correct transformation properties. We look at the example shown in 

Fig.n 




FIG. 1: A typical ampitude 



The amplitude is represented by an expression of the form: 

n"'(fcl) - j dk2 T''''{ki,k2,-ki-k2)D,d{~k2)T'"'f{~ki,~k2,ki+k2)D,f{ki+k2) (14) 

Using H12|l and (|13|l it is easy to see that this self energy transforms according to: 

n-''(fci) = v'^,{k,)v''-^{-k,)w\k,) . (15) 

In the sections below, indices in the 1-2 basis will be denoted bi and take the values 1 or 2. Keldysh indices will be 
written and are assigned the values a = 1 := r and a = 2 := a. R/A indices will be denoted Xi and assigned the 
values X = 1 := R and X = 2 := A. An n-point function in the 1-2 basis will be written: 

G{ki,k2,k^, - ■ ■kn)bi,b2M,---b^ '■— Gbi,b2M.---b^ (16) 

An n-point function in the Keldysh basis will be written: 

An n-point function in the R/A basis will be written: 

G(fci, A:2, fcs, • ■ ■ fcn)xi,X2,X3,---x„ Gxi,X2,X3, - x,^ (18) 

Note that the conservation of momentum gives fc„ — — (fci -I- ^2 + • • • + fcn-i) so that each n-point function depends 
on only n ~ 1 independent momenta. 



IV. THERMAL FUNCTIONS 



We write the general distribution fmiction as /(fc) and define the symmetric and anti-symmetric combinations: 

TA{k)^ f{k)- f{-k)- Ts{k) = f{k) + f{-k) (19) 
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In the Keldysh representation it will be useful to define the functions: 

^{h; ai) = Sai,a^A{ki) - 5on,r^s{h) ^{kf, "i) = 6ai,r^A{ki) - 5ai,a^s{ki). (20) 

To illustrate these definitions consider the expression: 

T{ki\a-i)T{k2\ 02)^(^3; as)G{ki,k2, fc3)ai,c<2,a3 (21) 
If we take the case a\ = a2 = r; 0:3 = a we obtain: 

:Fs{ki)J^s{k2):FA{k3)G{kuk2,ks)rra- (22) 
For bosons in equilibrium these expressions become: 

m ^ n{k) = l/(e'5'=° - 1) (23) 
:F{ki; ai) N{ki; a,) ; :F{ki; ai) N{ki; ai) 
J'sih)^-!; J^Aiki) ^ N{ki) 



where we have defined: 



N{ki; ai) = 5a.,aN{ki) + S^^r (24) 
N{ki;ai) = 5a,.rN{ki) + J^^.a- 
Ni := N{ki) = 1 + 2n{ki) 



Equilibrium functions satisfy: 



Ni+N2^Q if = -ki (25) 
1 + TViTVa + A^iiVa + N2N3 = if fca = -fci - 

Ni+N2 + N3 + Ni + N1N2N3 + N2N3N4 + NsN^Ni + NiN3N2 = if A;4 = -ki - k2 - ks 



where the dots indicate that higher order expressions can be generated by iteration. (Expressions with the same 
properties can be defined for the fermion distribution function: n/(fc) = l/{e^''° + 1) ; Np{ki) = 1 — 2n/(fci)). 
In order to write the KMS conditions in a concise way we define the function 

n 

C„:=C({7VO) = (7(Ari,7V2,...iV„)=^-[l-(-l)"+f] ■ N{n;p) (26) 

p=0 

where the symbol N(n;p) means the following: 

[a] start with n indices {xi, X2,X3 • • • x„} 

[b] consider all possible subsets of these indices containing p < n of the Xi 's (without considering order) 

[c] for each of these subsets, take the product of the corresponding iV(a;i)'s 

[d] sum over all sets 

In addition, we define N(n; 0) = 1. A few examples will illustrate this notation. 

Example [1]: if n = 3 and p = 2 then the possible sets of p arc: {xi,X2}, {x2,X3}, {x3,xi} and the result is 
N(3; 2) = N1N2 + N2N3 + N3N1. 

Example [2]: if n = 4 and p = 3 then the possible sets of p are: {a;i, a;2, 3:3}, {x2,X3,X4}, {a;3,a;4,a;i}, {x4,xi,X2} 
and the result is N(4; 3) = N1N2N3 + N2N3N4 + N3N4N1 + N4N1N2. 

Example [3]: if n = 3 and p = 1 then the possible sets of p are: {xi}, {X2}, {X3} and the result is 
N{3;l) = Ni + N2 + N3. 

Below we write out the first few C{Ni,N2, ■ ■ ■ -A^n)'s: 

C{m) = 1 (27) 

C(iVi,A^2) = iVi+A^2 

C{Ni,N2, N3) = l + N1N2 + N2N3 + N3N1 

C{Ni, N2, N3, N4) = Ni + N2+N3+N4 + N1N2N3 + N2N3N4 + N3N4N1 + N4N3N2 
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Note that because of (|25|l each of these expressions is zero if the momenta satisfy ki + k2 + ■ ■ ■ kn = 0- Using (|24|l and 
(|25|1 we have: 

C{N„ . . . 7V„) = • (28) 

n [Ki + ■ ■ ■ + Kn) 

To write the KMS conditions for a given n-point function in a compact form, we will need to use C-functions of the 
form defined above, but with arguments that are not the full set {Ni} with i running from 1 to n. We will need the C 
variables whose arguments are a subset of the Ni's. When working in the Keldysh basis we want C-functions whose 
arguments are the set of Ni^s whose corresponding cti's take the value [i] r or; [ii] a. When working in the R/A ba- 
sis we want C-functions whose arguments are the set of iV^'s whose the corresponding Xi's take the value [i] R or; [ii] A. 

We define these modified C-functions as follows. In the Keldysh basis the set {cti} contains n variables. The number 
of r's is defined to be rir and the number of a's is defined to be n^. Of course, we have n = Ur + Ha- We construct a 
set of rir variables: {NiSa-^r} and a set of Ua variables: {NiSai,a}- Using these sets of TV's wc define the corresponding 
C-functions: 

C„, := Ci{N,d^^^r}) (29) 

When working in the R/A basis we make the analogous definitions: The set {Xi} contains n variables. The number 
of i?'s is nji and the number of v4's is ua- We define: 

C„, := C{{N,dx,,R}) (30) 
C„^ :=C(Wx„a}). 

To clarify this notation, wc look at the following example: Take n — 7 and consider the set {ai} = {r,r,a,r,a,a,r}. 
We have: Ur = 4; C„,. = C(A^i, iVa, iV4, N7) and Ua = 3; C„„ = C{N3, N5, Nq). 

V. THE KELDYSH BASIS 

The rotation from the 1-2 representation to the Keldysh representation is accomplished by using the transformation 
matrix: 



UKeldysh^(l-2) ^ ^ I -1 ) ' ("^"^^ 



Note that in this case (|12() gives V — U. In general the n-point function in the Keldysh representation is given by: 

Gai---Q„ = 2~^1 Uai^ ■ ■ ■ UaJ'"Gbi--b„ ■ (32) 

The factor 2"/^^'-'^ is introduced to produce simpler expressions for the vertices, following and 1^. For example, 
using the notation defined in H32|) . the bare 3-point vertex in the Keldysh representation is written in tensor form: 

{{0,-i},{-«,0},{-z,0},{0,-i}} (33) 

This expression can be compared with the bare 3-point vertex in the 1-2 representation: 

{{-z,0},{0,0},{0,0},{0,z}} (34) 

In addition, the appropriate components of the Keldysh vertex correspond directly to the retarded vertices (see for 
example Eqn. (|41|l 'l. 

A. 2-point Function 

The 2-point function in the Keldysh basis is obtained from H32|) : 
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and can be represented (using Q and Hill) ) by: 

J^Keidysh := ( n'^'' ) (^^) 

'^^^^ ' ~ V -^11 - £'21 yl ^ Dadv 

Note that Z?aa — 0. In fact, for any n-point function in the Keldysh basis, it is always true that 

Gaaa-a = . (37) 

To understand this point, recall that according to Q), © and ^ it is always possible to express one component of 
the tensor 2-, 3- and 4-point functions in terms of the other components. It is straightforward to show that in the 
1-2 basis a constraint of the same form as (@J), © and ||HJ) exists for any arbitrary n— point function. In the Keldysh 
representation this constraint takes the form l|37() . The vertex function obeys an analogous constraint: 

prrr-r^Q^ (38) 



B. 3-point Function 

The 3-point function in the Keldysh representation is obtained from (|32|l : 

Gaia2a3 = ^U^^'-U^^^U^^^ Gbib2b3 ■ (39) 

We give three examples of the results that we obtain by summing over indices and using ®: 

Grrr — Gm + 6*221 + G212 + G122 (40) 
Grra — G\\\ — G\Yi + G22I — G222 
Graa — G\\\ — G\\2 — G\2\ + G122 • 

Using jSJ it is easy to show in co-ordinate space that 

Graa{xx,X2,X'i) = Gb.A^1^ ^2, X3) (41) 

where Gr-^ is the component of the 3-point function that is retarded with respect to the first leg and is given by the 
familiar expression: 

Gb,^{xi,X2,X3) (42) 
= i-if (eih - t2)e{t2 - i3)([[<^(a;i), 0(a;2)], ^(Xg)]) + e{h ~ h)e{h - i2)([[0(xi), 0(X3)], 0(X2)])) 

= {-if 0{h,ta)e{ta.h){mh).<l^{ta)i4>{h)]). 

{ta,tb} = {t2,t3} 

The summation indicates that we sum over the combinations: {ta = t2', 1^=1^) and [ta = t^] tb = ^2)- The other 
two components given in (|40|l correspond to similar expressions in co-ordinate space, but involve anti-commutators. 

C. 4-point function 

The 4-point function in the Keldysh representation is obtained from (|32|l : 

We give one example of the results that we obtain by summing over indices and using ((SJ: 

Graaa = Gnu ~ G1112 — G1121 — G1211 + G1122 + G1212 + G122I — Gi222 (44) 
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In co-ordinate space this expression has the form: 

Graaa{xi,X2,X3,X4) = G Ri{xi , X2, X3, X4) (45) 

= i-if ^(ii,ta)e(ia, 4)^(4, ic)([P(ii),0(ia)],0(ifc)],0(ic)]> 

{ta,*6,tc} = {t2,t3,t4} 

where the sum is over all possible assignments of the variables {t2, t^, to the variables {ta, t^, t^} and the notation 
Gri{xi, X2, X3, X4) indicates the 4-point function that is retarded with respect to the first leg. Expressions for higher 
n-point functions are obtained similarly. 



VI. THE (R/A) BASIS 



The matrix that performs rotations from the 1-2 basis to the R/A basis is: 



(fl.M)^(l-2) 



-fik) fik) 
1 f{k)/f{-k) 



(46) 



We will find it more useful to obtain expressions in the RA basis by rotating from the Keldysh basis. The matrix that 
transforms from the Keldysh basis to the 1-2 basis is given by the inverse of H31() . Combining with (|46() we obtain the 
matrix that transforms from the Keldysh basis to the R/A basis: 



U(R/A)^Keldysh — U {^pi/ A)^{l-2){U K eldysh^(l-2)) ^ 



(47) 



~V2f{k) \ 



The n-point function in R/A basis is given by: 

1^'^Gx,-x^ = U(kx)^l ■ • •{/(fc„)£"G„,...a„ . (48) 
General expressions that relate R/A n-point functions and vertex functions to Keldysh functions are given below: 
G(fci,---fc„)xi...x„ (49) 
_ (-l)"n^[^x„A + ^x„fl2/(fc,)] 



2"-^ nj'5x..fl + 5x„A/(-A;. 

r(fci,---fc„)^i-^" 



ai^{r,a) \l=\ J 



iR + SxiA^air ^{kl^ ' ' ' k^)"^ 



From (I27I), (|2H1) and 1(33 we have 



Gr...r^Q; r^ -^ = o 



(50) 



These equations can be written in a simpler form in equilibrium. There arc two kinds of simplifications that occur. 
We can replace general distribution functions with thermal distribution functions and use identities of the form H23|l , 
(I24|l . H25|) . (|28|l for combining groups of distributions functions. We can also use KMS conditions. The KMS conditions 
are a set of equations which are valid only at equilibrium that relate various components of an n-point function or 
vertex function. These conditions will be derived in the next section. For completeness, we give below two general 
expressions: the first corresponds to H49|) using thermal distribution functons, and the second comes from H49I) using 
thermal distribution functions and the KMS conditions. The general equation obtained from (|49l) using thermal 
distribution functions is: 



2nH-l(_i)nH + l CnA{N^Sx,.A})Gik,,■■■kn)x,■■■X^ 

= Gn^{{NiSxiM}) ■ ^ I J]^ [iV(fc;;a;)(5x,A +^X,_R'5a,a ] G{ki, ■ ■ ■ kn)ai . a„ 
aie{r,a} = l 

2"^-^(-l)"-^+iG„,({iV,5x,,fl})r(fci,...fc„)^--^" 

,R + SxiA r(fci, ■■■k 

n ) 



(51) 
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After applying the KMS conditions to the right hand side of this expression we obtain: 

2n«-l(_l)n. + l G(fci, . . . fc„)x,...x„ (52) 

= ^ I JJ^ aO'^^ifl + (^^iA^Qia 1 G*(fci, • • • fc„)ai - Q„. 

aie{i\a} = l I 

For clarity we give a few examples. We first give general expressions using 149|l . We next give the results obtained with 
thermal distribution functions using 1)51(1 . Finally, we give the expressions obtained by applying the KMS conditions 
to the right hand side of the previous equations, or equivalently from using l|52|) . 

[i] The 2-point functions in the R/A formalism are given by: 

^RR = (53) 

DrA -^s{k2)Dar 
DaR = -Ts{ki)Dra 

DaA = {Ts{ki)Ts{ki)Drr - TA{kl)Ts{k2)Dar - Ts{ki)TR{k2)Dra) I (2/(-fci)/(-fc2)) . 

Using equilibrium distribution functions these results reduce to: 

Drr^Q] DRA^Dar] DAR^Dra] D AA ^ 2n(-fci)n(-fc2) '^^'•'•^^^ + + ^■^^r<^) ■ 

Applying KMS conditions we have: 

Drr = ; Dra -D*a ; Dar. = D*^^ ; Daa = , 
which can be written in the familiar form: 

( Drr Dra \ f Dadv \ /r 

°«/^ = [ Dar Daa,) = [ Dr., ) ' (^4) 

[ii] The 3-point functions in the R/A formalism are given by: 

GrrA = T!S4^^s{ks) Gaar (55) 

f(ki + k2) 

^ _ /(fcl)/(fc3) ^ I, . ^ 

f{ki + ka) 

^ /(fc2)/(fc3) ^ . ^ 

J[k2 + kz) 

Gaar = -\ / ^ {J^s{kl)J^s{k2)Grra - J^A{ki)J^s{k2)Gara - ^siki)TA{k2)Graa) 

2 /(-fcl)/(-K2) 

GaRA = ~\ {J^s{kl)Ts{k'i)Grar - ^A{ki)Ts{kz)Gaar - ^ s{kl)T A{k:i)Graa) 

GrAA = /} 1^1 J-, {^s{k2)J^s{k'i)Garr - J^A{k2)Ts{kz)Gaar - -^S (fc2 )-?^A (fc3)Gara ) 

2 J(-k2).f[-k3) 
1 

Gaaa 



4 /(-fci) /(-fc2) /(-fcs) 

J^sikl)J^s{k2)^s{k3)Grrr 
-TA{kl)Ts{k2)Ts{k3)Garr - (^1 )-?^A (^2 (^3 ) G^ar " (^1 (fe )-^A (^s) G^ra 
+TA{kl)TA{k2)TA{k3)Gaar + A{kl)T s{k2)T A{kz)G ara + J^s{kl)TA{k2)TA{kz)Graa 

Using equilibrium distribution functions these results reduce to: 

GrrA = —- ^—Z ^ Gaar (56) 
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Grar 
Garr 
Gaar 
Gara 
Graa 
Gaaa 



(iVi + iVg) 



{N2 + N3) 



G„ 



1 
1 

{N2 + N3) 



{Grra + N\Gara + ^2Graa) 
{Grar + N\Gaar + N^Graa) 



{Garr + ^2Gaar + ^sGara) 



1 



4 n(-fci) n{-k2) n^-k^) 
Applying KMS conditions we have: 



{Grrr + N^G^rr + A^2Grar + N^Grra + NiN^Gaar + NiNsGara + N^NsGraa) 



Grra 

Grar 

Garr 

Gaar 
Gara 
Graa 
Gaaa 



2 
1 

~2 
1 

~2 

aa 



^raa 



{G*arr + ^^G*^^^ + N^G^^^) 



= 0. 



For the 4-point functions using equilibrium distribution functions we have: 
C{Nr,N2,N:,) 



Grrra 
Grrar 

Grarr 
Garrr 

GrraA 

Grara 
Graar 
Garra 
Garar 
Gaarr 
Gaaar 
Gaara 
Garaa 



C{NuN2,Ni) 
4 

C{NuN3,Ni) 
4 

C{N2,N3,N4) 

4 

l C{Ni,N2) 
2C{m,Ni) 
l C{NuN3) 
2C{N2,Ni) 
l C{Ni,N4) 
2C{N2,N3) 
1 C{N2,N3) 
2C{NuN^) 
1 C{N2,N^) 
2C{Ni,N3) 
l C{N3,Ni) 
2C{NuN2) 
1 

C{NuN2,N3) 
1 

C{NuN2,N4) 
1 

C{NuN3,N4) 



(57) 



Gaaar 
Gaara 
Garaa 
Graaa 

{Gaarr~^^zGaaar~^^AGaara) 

{Garar~\-^2Gaaar-\-^4Garaa) 

{Garra-^^2Gaara-^^zGaraa) 

{G raar'^^lG aaar"^-^ aG raaa) 

{Grara~^^\ Gaara~^^Z Graaa ) 

{Grraa-\-Nl Garaa-^^2 Graaa ) 
{Grrra~\~ ^lGarra~\~ ^2G 

rara + NsGrraa+mN2Gaara+NiNsGaraa+N2N3Graaa) 

(G 

rrar 

+ NiGarar+N2G 

raar rraa 

+ N,N2Gaaar+NiN4G 

araa + N2NiGraaa) 

( Grarr+ N\ Gaarr+ N'iGraar+ NiGrara+ ^1 -^3 Gaaar+ -^1 -^4 Gaara+ N^NiGraaa ) 



(58) 
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GrAAA = J.J TZ TfT- [G arrr+ N2G aarr+ N^G arar+ NiG arra+ N2NzG aaar+ N2N4^G aara+ N2,N iO araa) 

G(iV2, -/V3, iV4j 



^^^^ g n(— fci) n{~k2) n{—k^) n{—k/i) 
-I-/V, /V„r7 -I- /V, /V„r7 4^ N, 1 



Applying KMS conditions we 



+ N2G rarr + N3G rrar rrra 

+N^N2Gaarr + N^N^Garar + N.N^Garra + A^2A^3G™ar + N2N4Gra^ 
+N^N2NsGaaar + N^N2NiGaara + N^N^N^Garaa + N2N^NiGraaa) 

have: 



rara 



■ N^N^Grraa 



Grbra 


^ ^G ^J^^^Q^ 


AT ^* 1 AT ^* 1 7y /^* 




GfiJfAB, 




1 ^ arar^ raar^ ^ rraa 


- /Vi /Vo(7* + /Vi A^^ a* + /Vo A^^ (7* 


G BARB, 




\^ aarr^ 'A raar^ ^ rara 


-NtN-^G* -^NtNaH* -^N-^NaH* 


Garbb. 




JV C* -\- JV C* 4- N C* 

aarr^ arar'^ ^ arra 


TV T TV T 1 TV T" AT" i TV T AT" 

'N2N3G^^^^+N2N4G^^.,^+N3N4G^^^^ 


GrraA 






Grara 


= ^2 (^rara+^lGaara+^sG*^,^,^) 




'^RAAR 


= ^2 (^™<"-+^l^aaar+^4G*^^„) 




Garra 


= -i (G* +A^9G* +7VciG* ) 
2 V^arra' -^^ aara ' 0'-' araa/ 




Garar 


= -1 fG* +N9G* +NaG* ) 

2 arar ' ' z'-' aaar ' 'J:^ araa/ 




Gaarr. 


V^aarr' o'-' aaar 1 ' 4'-' aara J 




Gaaar. 




^ aaar 






Gaara 




^ aara 






Garaa 


_ 






Graaa 








Gaaaa 


= 







(59) 



VII. KMS CONDITIONS 



The KMS conditions are a set of relations between the various components of a given n-point function that hold 
at equilibrium. The KMS conditions are often useful for simplifying the expressions that result from the contractions 
over indices, in any representation of real time statistical field theory. The KMS conditions have the simplest structure 
when expressed in the R/A basis. 



A. Keldysh Basis 

The KMS conditions in the Keldysh representation can be written: 

G„„({iV,(5„;,a}) (f[[^(M'^(^l)^c.lr+Sc.'^aSa,a 
aie{r,a} \l = l 

= Cn^{{N,6a,^r}) (f[[N{ki;ai)6^>a+Sa[rSa,a 

ai<£{r,a} = l 



Gax-'-CXn 



n 



(60) 



In this master equation, the set of variables {a'^} are external variables that are not summed over. Thus, in principle, 
(|6(JII contains 2" equations which come from the 2" choices of the sets {a'^}. In fact, one half of these equations is 
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the complex conjugate of the other half and thus we have 2" ^ KMS conditions for each n-point function. We give 
several examples below: 



[n = 2] 

For the 2-point function we obtain 2"~^ 



n=2 



= 2 equations: 



[n = 3] 

There are 2"-i 



n=3 



Dra {P) = Dl, (p) or Dret ip) = Dl^, (p) 
Drrip) = N{p){Dra{p) - £>ar(p)) Or Dsymip) = N{p){Dret{p) - Dadv{p)) ■ 

4 independent KMS conditions which are: 



(61) 



Grrr + N^Garr + N2Grar + N^Grra + NiN2Gaar + NiN^Gara + N2NsGraa = 

Grra + ^l^a™ + N^Graa = {N, + N2)Gl,r 
Grar + N^Gaar + N^iGraa = {Nl + N^)Gl,.^ 
Garr + N^G^ar + N^^Gara = (^2 + N^)G;^a 



(62) 



[n = 4] 

There are 2"-^ 



n=4 



= 8 independent KMS conditions which are: 



rrrr arrr rarr + N3G rrar rrra + N1N2G a 

+ N2N3Graar + N2NiGrara + N^NiGrraa + NiN2N3Gaaar 
+ NiNaN^Garaa + N2N3NiGraaa = 



■r + NiNsGa 
NiN2NiGa. 



N^NiGa 



(63) 



Grrra+NiGarra+N2Grara+N3Grraa+NiN2Gaara-\-NiN3Garaar\-N2N3Graaa 

= (1 + iViiVa + NiNs + ^2^3)G:„„, 

Grrar'^ Garar + N2Graar+NiGrraa+NiN2 Gaaar+ -^1 NiGaraa+ N2N4Graaa 
= (1 + N1N2 + NiNi + N2Ni)G:,,, 

Grarr+ NiGaarr+ N3Graar+ NiGrara^ NiNsGaaar+ NlNiGaarar\- ^zNiGraaa 

= (1 + N.N:, + NiN^ + iV3^4)G:,„„ 

Garrr+ N2G aarr+ N3G arar+ NiGarra^ N2NsGaaar+ ■^2-^4Gaora+ N^NiGaraa 

= (1 + iVaiVs + N2Ni + N3N^)Gl,,, 



{N3+N4) (G 

rraa araa 

+N2G 

r 

{N2+N4) (G 

rara aara 

+N3G 

r 

{m+N4) (G 

arra 

+N2G 

aara a 



,) = {N,+N2) {G:,rr+NsG:,,,+N,G:,,,) 

^) = {N1+N3) {Glrar+N2Gl,,,+NiGl,^,) 

^) = {N2+N3) {g;,,,+NiGI,,,+N4g;,^j 



[n = 5] 

There are 2"-i 



16 independent KMS conditions. We write down three representative ones: 

N.Gr 



n— 5 

Grrrrr N\G arrrr H~ ■N2Gj-arrr H" -^zGrrarr N^.Gfrrar ^ -^^Grrrra 

+ N1N2G 

aarrr 

+ NiN^Gararr + N.N^G 

arrar 

+ NiN^Garrra + N2N3G. 

rarar + N2N5G 

rarra + N3NiG rraar + N3N5G 

rrara 

+ NiN2N3Gaaarr + NiN2NiGaarar + iViTVaiVsGaarra + NiNsNiGaraar 
+ NiNsN^Garara + NiN^N^Garraa + N2N3NiGraaar + N2N3N5Graara 
+ N2NiN5Graraa + NzN^N^Grraaa + NiN2N3NiGaaaar + N1N2N3N5G. 
+ NiN2NiN5Gaaraa + iViiV3iV4 A^gGa^aa + A^2iV3iV4iV5G™aaa = 



(64) 



raarr 
rrraa 



Gf 



N,Ga 



iVoG, 



N3Gr 



N^Gr 



NiN2Ga 
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+ NiN^Garara + NiN^Garraa + N2N'g,Graara + N2N4,Graraa + N^N^Grraaa 

+ NiN2N3Gaaara + NiN2NiGaaraa + NiN^N^Garaaa + N2N3.NiGraaaa 

= {Ni +N2+N3+Ni+ N1N2N3 + N1N2N4 + NiNsNi + 7V27V3^4)GLaar 

(A'4 + A^5) {Grrraa + NiGarraa + N2Graraa + N^Grraaa + NiN2Gaaraa + NiN^Gar 

= (1 + N,N2 + N.Ns + N2N3) (GL,„, + N^G:,,,, + A^sGL^™) 



N2N3G. 



raaaa ) 



The KMS conditions for the vertex functions have almost exactly the same form. The master equation is obtained 
from (|60|) by interchanging the indices r and a: 



a; £ {r.a} \/ — 1 



r^air 

71 



(65) 



B. R/A Basis 

The KMS conditions have an even simpler form in the R/A representation. This can be anticipated by comparing 
the general structure of (|51|l and 1601) . The general expression for the KMS conditions in the R/A basis is: 



(66) 



This master equation contains a series of equations where the variables Xi take all possible combinations of the 
values {R, A}. We use the notation: A = R, R = A so that = ua and = nn. We give several examples below: 



= 2] For the 2-point functions we have: 



DABip)=D*j,Aip)-, Daa = 0- 



= 3] For the 3-point functions we have: 



Gaaa — 
2CiN3)GRRA = - G(7Vi, N2) G\j,^ 
2C{N2)Grar. = - CiNi, N3) G*ARA 
2C{N^)Garr = - C{N2, 7V3) G^AA 



(67) 
(68) 



'^GjfjiA 
2GiiAR 

2Gab.r 



{N,+N2)G*^^j, 
{N,+N3)G*ARA 

{N2 + N,)G*RAA 



where the arrows indicate the results that are obtained by replacing the G-functions by their definitions in terms of 
thermal functions (see Eqn. (|26|l ). 

[n — 4] For the 4-point functions we have: 



Gaaaa = 
CiN3,Ni)GRRAA = C{N,,N2) G\arr 
C{N2,Ni)GRARA - G(7Vi,iV3) G\j,Aj^ 
C{N2,N3)Graar = G(A^i,iV4) G*ARRA 

^GaRRR = C{N2,N3, Ni) G*RAAA 
4GrarR = G(7Vi, 7V3, Ni) G*ARAA 

4Grrar = GiNi,N2, Ni) G\j,j,j, 
^Grrra = GiNi,N2,N3) G\j,j,j, 



(69) 



(iVa + Ni)GRRAA = {Ni + N2) G*AARR 

{N2 + Ni)GRARA = {Ni + N3) G\j^AR 

{N2 + Nz)Graar = {Ni + Ni) G*ARRA 
^Garrr = (1 + N2N3 + NsNi + NiN2) G*raaa 

4.GRARR = (1 + + N^Ni + N^Ni) G*ARAA 

'iGRRAR - (1 + 7Vi7V2 + iV2A^4 + N^m) G\ara 

AGrrra = (1 + A^1^2 + N2N3 + N3N1) G\aar 



5] For the 5-point functions we give a few examples: 



Gaaaaa — 
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C{N3,m,N5)GRRAAA 

^ (1 + NsNi + NiN5 + N5N3)Grraaa 

8GRRmiA 

— * ^Grrrra 

Note that in every case one of the KMS equations requires the vanishing of the Ga - a component of the n-point 
function. The rest of the KMS equations are simple relations between pairs of off-diagonal components. 

The KMS conditions for the vertex functions have almost exactly the same form. They are obtained from the same 
master equation, with the indices R and A interchanged: 

2"--"- c„j{iv,<5x..K}) r^^ -^" = {-iTc^nm5^.,R}) (r )^- -^" (7i) 

VIII. A PROGRAM TO PERFORM THE CONTRACTION OF INDICES 
A. Description of the Program 

Our program calculates the integrand for any diagram by contracting indices. The calculation can be done in the 
1-2, Keldysh or R/A basis. The basic strategy of the program is to treat a Feynman diagram as a tensor product of 
propagators, and bare and corrected vertices. In principle, this type of calculation can be done by hand, but when 
working with more than a few indices the process becomes extremely tedious. The symbolic manipulation program 
Mathematica is ideally suited to perform this kind of tensor calculation. 

The program is divided into six main sections. Only the Input section needs to be edited by the user. The user 
enters some variables in this section, in order to specify the diagram that he wants to calculate. The rest of the 
program can be immediately executed, and the result is output in the last section. The user has the option to output 
the results to a file. The results are also defined functionally. Some basic functions are defined within the program 
that can be used interactively to manipulate the result. 

When working in the RA basis the program always assumes thermal distribution functions, and thus the result is 
only valid in equilibrium. The user has the option to implement the KMS conditions, since both forms of the result 
can be useful. In the Keldysh basis, thermal distribution forms are only assumed if the KMS conditions are used. In 
the 1-2 basis the KMS conditions cannot be implemented in the current version of the program. Note that we have 
not used different notation to indicate equilibrium and non-equilibrium distributions (see Eqns. I|23|) and 124|) '): the 
output of the program will always contain distribution functions written in the form Hp or Np. 

The program is designed to contract CTP indices and is written for scalar bosons. We set all coupling constants to 
one and use the notation: 





= ip 


> 




X 






= -in 


> 


= -ir 


X 


= -iM 


:^ 


= -iC 



FIG. 2: Definitions of notation for propagator and vertices 

Note that additional numerical factors are introduced in the Keldysh representation, as defined in (|32() . The program 
can be used for other field theories with appropriate adjustments. The appropriate coupling constant (s) must be 
inserted by the user. In addition, any dirac, lorentz, or other group structure must be separately handled by the user. 

We describe briefly the main sections of the program: 
1. The Initialization section inputs the necessary Mathematica package. 



= ~2G{N,,N2)G*^^j,^pj, 

= -2(7Vi + N2) G\a„p„ 

= -C{N,,N2,N3,N^)G*aaaab^ 

= -(iVi +N2+N3+N4, + N1N2N3 + N2N3N4, + N^N^Ni + N4N3N2) G*aaaar 
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2. The Input section is edited by the user to input specific parameters corresponding to the diagram he wants to 
calculate. This process is described in detail below. 

3. The Definitions section establishes some basic definitions that will be used throughout the program. 

4. The Find Loops section identifies all closed loops in the diagram. The variable {loop zero subs) is printed out. 
The usefulness of this variable is explained in point 11 below. 

5. The Calculate Diagram section performs the calculation. 

6. The Results section outputs the results of the calculations and performs some basic manipulations to simplify 
them. 



B. Input section 

We illustrate the Input section of the program with an example. Consider the diagram in Fig. O 



Q 




FIG. 3: Example Feynman diagram. Uppercase letters are momenta and lowercase letters label indices of the propagators and 
the vertices. 

The Input section for this diagram is reproduced on the page that follows the itemized list below. Each entry is 
explained in the order it appears in the Input section. 

1. number of external legs, external momenta, external indices: Specify the number of external legs, 
and the momentum and corresponding index for each leg. 

2. number of internal indices, internal indices: Specify the number of internal indices and list them. 

3. list of momenta, list of indices: For each bare vertex, list the incoming momentum for each leg and the 
corresponding index. 

4. list of momenta, list of indices: For each corrected 3-point vertex, corrected 4-point vertex and corrected 
5-point vertex, list the incoming momentum for each leg and the corresponding index. 

5. list of momenta, list of indices: List the momentum and the corresponding pair of initial and final 
indices for each propagator. 

• In this example, there is one corrected 3-point vertex and one corrected 4-point vertex. The corresponding 
momentum arguments are listed in double set brackets, indicating that the one vertex is the first in a list of 
length one. The same notation is used in any case where multiple entries would create nested lists. This includes 
the momenta and indices of bare or corrected n-point vertices, with n > 3, and the indices of propagators. The 
momenta for any number of propagators (even if there is only one) are listed in single set brackets. 

6. choose basis (onetwo, Keldysh or RA): Indicate the basis as either onetwo, Keldysh or RA. 

7. combination(s) to be evaluated: ie. ra, ... or AR, . or All: Indicate the set of external indices 
of the diagrams that are to be calculated or simply specify All. 
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8. Simplify result? (option to use the Mathematica function 'Simplify'): The user can choose to 
have the program apply the Mathematica function Simplify to the result [simplifylt = yes). For smaller di- 
agrams this produces neater results. For larger diagrams it can lead to significantly longer running times, 
without producing a result that is much more compact. 

9. use the KMS conditions: In the Keldysh or RA bases, one can choose to enforce the KMS conditions [useKM- 
S conditions = yes). In the 1-2 basis choosing {useKMSconditions=yes) has no effect on the output. If the user 
has chosen the Keldysh basis, he can enforce the KMS conditions for the propagator, but not for higher n-point 
vertices {removeFs=yes) . Note that if the user chooses useKMSconditions=yes and removcFs—no, the first 
choice over-rides the second, and KMS conditions are used for all n-point functions including propagators. 

10. express answer in terms of RA expressions: In the Keldysh basis, if the calculation was done using the 
KMS conditions, the user can request that the right-hand side of the final expression be written in the R/A 
basis {InTermsOfRA = yes), since this basis frequently produces more compact results. The default value of 
this parameter is "no" and produces results in which both sides of the equation are written in the same basis. 

11. remove terms which are zero after integration: One can choose to eliminate terms which will vanish 
after integration (removeZeros—yes). This option only works in the Keldysh or RA bases (in the 1-2 basis 
choosing (removeZeros—yes) has no effect on the output). The program identifies closed loops in the diagram, 
and then looks for terms in which all of the propagators that form a closed loop have poles on the same side of 
the real axis. Choosing [removeZeros^yes) will cause the program to set all of these terms (except for tadpoles) 
to zero. The program automatically exempts tadpole loops, so that they are not incorrectly set to zero. The list 
of combinations of propagators which will be removed by this option are stored in (loopzerosubs) and printed in 
the Find Loops section of the program. The propagators that correspond to tadpoles are also listed. One can 
compare these lists to the original diagram as a check of the information given in the Input section. For our 
example we have: 

G,4L)Gr4L + P + Q)^0; G^,{L)G^,{L + P + Q) ^ Q (72) 

The two expressions correspond to the one loop in our example diagram, with the momentum routed clockwise 
or counter-clockwise. 

12. replace C functions in terms of N functions: If the KMS conditions are used, the program does the 
calculation in terms of the C-functions defined in Eqn. (|26(l . One can choose to replace these C-functions with 
their definitions in terms of thermal fmictions [replaceCs^yes) - see Eqn. H26|l . 

13. specify nsune of output file using quotes: The user can specify a filename to output the results to, or 
specify an empty string ("") for no output file. 

14. Directory for output file is directory of this notebook as specified below: The directory of the 
output file is set to be the same as the directory of the original notebook. The user can specify another directory 
for the output file by rewriting this line. 



C. Error messages 

Some checks have been implemented to detect possible input mistakes, and give error messages indicating the nature 
of the problem. A few examples are given below: 

1. If the incoming momenta for any of the internal vertices, or for the external legs, do not sum to zero, the error 
message appears: "Momentum non-conserving vertex detected." 

2. If the combinations of external indices requested using "combination(s) to be evaluated" do not correspond 
to the basis specified using "choose basis (onetwo, Keldysh or RA)," an error message appears. For ex- 
ample, choosing the RA basis and requesting the combination rra (corresponding to the vertex function P''"'^) 
produces the error message: "Initialization Failed: Combination rra is not specified in the RA basis." 

3. If the combinations of external indices requested do not have the same number of variables as the number 
of external legs, an error message will appear. For example, calculating a 3-point vertex function (as in our 
example) and requesting the combination rraa (corresponding to the 4-point vertex F''''"") produces the error 
message: "Initialization Failed: Combination rraa does not have the correct number of indices." 
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Input 



Off [General : : spell] ; 
Off [General: :spelll] ; 

(* number of external legs: *) 

(* external momenta: *) 

(* external Indices: *) 

(* number of internal indices: *) 

(* internal indices: *) 

(* list of momenta: *) 
(* list of indices: *) 

(* list of momenta: *) 
(* list of indices: *) 

(* list of momenta: *) 
(* list of indices: *) 

(* list of momenta: *) 
(* list of indices: *) 

(* list of momenta: *) 
(* list of indices: *) 

(* list of momenta: *) 
(* list of indices: *) 



nvmiExternalLegs = 3; 
externalMomenta = {P, Q, -P-Q}; 
externallndices = {a, b, c}; 
numlnternallndices = 4; 
internallndices = {d, e, f, g}; 

bareSptMomenta = { } ; 
bareSptlndices = {}; 

bare4ptMomenta = { } ; 
bare4ptlndices = {} ; 

Cor3ptMomenta = { {P + Q + L, -L, -P-Q}}; 
CorSptlndices = {{f, g, c}}; 

Cor4ptMomenta = {{P, Q, -P - Q - L, L}}; 
Cor4ptIndices = {{a, b, d, e}}; 

CorSptMomenta = { } ; 
CorSptlndices = { } ; 

propagatorMomenta = {L, P + Q + L}; 
propagatorlndices = {{g, e}, {d, f}}; 



(* choose basis (onetwo, 

Keldysh or RA) *) 
(* combination (s) to be evaluated: ie. {ra, ...} or {AR, . } or 

All *) 

(* Simplify result? (using Mathematica function: 
Simplify) *) 

(* use the KMS conditions *) 
useKMSconditions = no; 

(* use KMS conditions for propagator only 

(in the Keldysh basis) *) 
(* express answer in terms of RA expressions 

(in equilibrium and Keldysh only) *) 

(* remove terms which are zero after integration *) 
(* replace C functions in terms of N 

functions *) 
(* specify name of file to output to using 

quotes *) 
FileName = "leafA_exainple.txt"; 

(* Directory for output file is directory of this notebook as specified below *) 
SetDirectory[ DirectoryName[ ToFileName[NotebookInformation[InputNotebook[] ] [[1, 2J ] ] ]; 



Basis = Keldysh; 
combinations =A11; 
simplifylt = no; 

removeFs = yes; 
InTermsOfRA = no; 
removeZeros = yes; 
replaceCs = yes; 
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4. Each external index must appear one time only as a vertex index. Each internal index must appear one time as 
a vertex index and one time as a propagator index. If the indices are incorrectly entered, an error message will 
appear. 

5. Because of the fact that vertices are defined with incoming momenta, for each propagator, the momentum 
corresponding to the second index must be the same as the momentum for the vertex leg with the matching 
index. Similarly, the momentum corresponding to the first index must be the negative of the momentum for 
the vertex leg with the matching index. If the momentum variables are incorrectly entered, an error message 
will appear. 

D. Interactive Manipulation of the Results 

The result of the calculation is stored in the internal variable diagram^ The user can extract and manipulate the 
results in several ways. 

1. One can extract a specific result using the function DIAGRAM. For our example, evaluating DI AGRAM[a^ r, r] 
gives the result for F"'"'" in the Keldysh representation: 

-UgMG^JL + P + Q)F-'^-(P, Q,L,^L^P- Q) (73) 
paar(^ + P + Q, -L, - P - Q) - N lT^'' {L + P + Q, -L, -P - Q) 
+ Nl+p+qT'^'{L + P + Q, -L, -P - Q) 

(the prefactor / j^Pyi explicitly written). 

2. In Eqn. (|77|l we define notation for the Keldysh basis so that each string of indices of the form rrar ■ ■ ■ 
corresponds to a single numerical index. A function called GAMMA is defined that will automatically supply 
the result of the calculation for a given numerical index. All vertices are translated into the notation defined in 
Eqn. lf77|) . Propagators are written using the notation Dra{P) — P>ret{P) '■= r{P); Dar{P) = Dadv{P) ■= a{P)- 
The function GAMMA is illustrated below. From Eqn. (|77fl . the index '2' corresponds to the string 'arr' for a 
3-point vertex function. For our example, evaluating GAMMA[2] produces the translation of the result for F"'''' 
that is given in ifTHI) : 

-h,a{L)r{L + 5)M(2, P,Q,L,-L- S) (74) 
F(4, L + S, -i, -S*) - A^lF(2, L + S, -L, -S) + Nl+sT{3, L + S, -i, -S*) 

IX. THE QED WARD IDENTITY FOR THE 4-POINT FUNCTION 

In this section, we present a calculation of the ward identity for the QED 3-point vertex function, and the QED 
4-point vertex function involving two photons and two fermions. Throughout this section we work in the Keldysh 
representation. As explained in the beginning of section rVIIII the goal of our program is provide compact results for 
the expressions that result from summing over CTP indices. The dirac structure of fermions and any group structure 
associated with all fields, must be handled separately by the user. For the calculations done in this paper, there are no 
additional complications. The ward identities are derived by comparing groups of integrands, without evaluating the 
integrals themselves. Consequently, we can simply suppress all dirac and lorentz indices. For example, the result for the 
contribution to 11(2, p) shown in Fig. (Sjis given 

bylHS)- The first term is -f / j^a{l)r{l+p)T{2,p,l+p)T{A,l+p,p). 
Including all indices, this factor gives a contribution: 

'I J (^«M-(0r^/5'(2,z,z + p)'^/3'a'(^ + p)r^'J4,; + p,p) (75) 

where {/x, j/} are Lorentz indices and {a, a' , P, P'} are Dirac indices. Throughout the rest of this section we will 
leave all integrands in the form described above, with dirac and lorentz indices suppressed. In addition, we define 
k + p + q = k + p = t, p + q = s and use the shorthand notation: Dret{p) = r{p), Dadv{p) = o,{p), Dsym{p) — f{p)- 
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In order to simplify the notation for the vertices, we replace each combination of the indices {r, a} by a single 
numerical index: 

. . .p„) ^ r(")(z,pi,p2, ■■■Pn) (76) 

We assign the choices of the variables 01102 • ■ ■ On to the variable i using the vector 

^"=(::)-«(2)«(:;) <"> 

where the symbol (8) indicates the outer product. For each n, the ith component of the vector corresponds to a list of 
variables that is assigned the number i. To simplify the notation we drop the subscripts and write a list like rir2a3 
as rra. For clarity, the results are listed below. 

[a] 2-point functions: rr — > 1, ar ^ 2, ra ^ 3, aa ^ 4 

[b] 3-point functions: rrr —^ 1, arr — > 2, rar — > 3, aar —^ 4, rra —^ 5, ara — > 6, raa 7, aaa —^ 8 

[c] 4-point functions: rrrr 1, arrr 2, rarr — > 3, aarr — > 4, rrar 5, arar — > 6, raar — > 7, aaar — > 8, rrra — s- 9, 
arra —f 10, rara — > 11, aara ^ 12, rraa 13, araa — ^ 14, raaa 15, aaaa 16 

[d] 5-point functions: rrrrr 1, arrrr 2, rarrr — )■ 3, aarrr 4, rrarr 5, ararr — > 6, raarr — 7, aaarr — 8, 
rrrar — > 9, arrar — > 10, rarar 11, aarar 12, rraar — )■ 13, araar — )■ 14, raaar 15, aaaar — > 16, rrrra — > 17, 
arrra ~> 18, rarra — + 19, aarra — > 20, rrara — > 21, arara — ^ 22, raara ~> 23, aaara — > 24, rrraa — > 25, arraa — ^ 26, 
raraa — > 27, aaraa — s- 28, rraaa 29, araaa — > 30, raaaa 31, aaaaa — > 32 

We note that i = 1 corresponds to a vertex that is identically zero, for any number of external legs, as a consequence 
of the general relation p'''''' ' '' = which can be obtained from (|65|l . 

To reduce the number of indices we will introduce separate names for the first five vertex functions and write: 

p(2) ^ J. . p(3) ^ p . p(4) ^ . p(5) ^ ^ (78-) 

The notation for the momentum arguments is as follows: 

S(Pj„) ; r>'{ptn,Pout) ; M'"'{p,n,qi,q2,Pout) ; C^"''^{p„i,qi,q2,q3,Pout) , (79) 

where pin is the momentum of the incoming fermion, {qi, q2, q^} are the momenta of the incoming photons, and pout 
is the momentum of the outgoing fermion. Note that the momentum of the photon is not written for the 3-point 
vertex since it can be inferred from the momenta of the fermions. Similarly, the self energy is written S(pi„) instead 

of E{pin,-Pin)- 

We begin by calculating the ward identity at the bare 1-loop level. In the next section, we verify that the same 
ward identities are satisfied by the complete set of graphs involving full corrected vertices. 

A. Bare 1-loop Diagrams 

We start by looking at 1-loop diagrams with bare vertices, in order to determine the form of the ward identities. 
For the 2- and 3-point functions the graphs are shown in Figs. (jS^) and ^^). For the 4-point function we have the 
box graph (Fig. (0;)), and the crossed version of the box graph where the two external photons are interchanged. For 
the five point function, the basic graph is show in Fig. ((SJi) • There are six versions of this graph which correspond to 
the six possible permutations of the three external photons. 
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FIG. 5: One loop diagrams 

We use the notation defined in Fig. |21 and insert the appropriate coupling constants for QED. We include the 
numerical factor as defined in H32|l : 

— zE^^-' ~ self — energy graph (80) 
— zeF^^-* — \/2 ^triangle graphj 

— ze^M*-^-* = 2 ^box graph + crossed box graph^ 

— ie^C^^-* = 2 \/2 ^six permutations of 5 — point graph^ 

where the superscripts refer to the loop order of the graph. The ward identities are of the form 

Q^r« (z, p, p + q) - (ji , p) - Ed) (j2 , P + g) (81) 
g'^M« {i,p,q,k,p + q + k)^ r(.i) in ,P,P + k)- rl^) (j2 , p + g, fc + p + g) 

Q^CW^ ii,p,q,k,s,p + q + k + s) = A41^ iji,p,k,s,k+p + s)- M W {j2,p + q, k, s, k + p + q + s) 

where the indices {i,ji,j2} refer to the choices of the variables aia2 ■ ■ ■ ctn as defined in (|77|l . We list below the sets 
of these indices for the 3-point, 4-point and 5-point functions. 

r: (2,2,2), (3,3,2), (4,4,1), (5,3,3), (6,4,4), (7,1,4), (8,2,3) (82) 

M: (2,2,2), (3,5,2), (4,6,1), (5,3,3), (6,4,4), (7,7,4), (8,8,3), (9,5,5), (10,6,6), (11,1,6), 
(12,2,5), (13,7,7), (14,8,8), (15,3,8), (16,4,7) 

C: (2,2,2), (3,9,2), (4,10,1), (5,3,3), (6,4,4), (7,11,4), (8,12,3), (9,5,5), (10,6,6), (11,13,6), 
(12,14,5), (13,7,7), (14,8,8), (15,15,8), (16,16,7), (17,9,9), (18,10,10), (19,1,10), (20,2,9), 
(21,11,11), (22,12,12), (23,3,12), (24,4,11), (25,13,13), (26,14,14), (27,5,14), (28,6,13), (29,15,15), 
(30,16,16), (31,7,16), (32,8,15) 

We give three specific examples below. The first set of numbers in the first line of H82|l is (2,2,2). Using (|77|l and H81|) 
the corresponding ward identity is: 

Q'^rW(2,p,p + g)-E(i)(2,p)-E(i)(2,p + (7) (83) 
or Q ■ r W {p, p + q) = E W {p) - E W {p + q) 

The first set of numbers in the second line of (I82|) is (2,2,2). The corresponding ward identity is 

Q^MW, {2,p,q,k,p + q + k) = rW (2, p, p + fc) - r W (2, p + g, fc, fc + p + q) (84) 

or Q^-ImI^I) (p,p + fc)=(rW) (p,fc+p)-(rW) {p + q,k+p + q) 

V / arrr V / arr V / arr 

The last set of numbers in the last line of H82|l is (32,8,15). The corresponding ward identity is 

,p,q,k,s,p + q + k + s) = Mj;,'^\8,p, k,s,k + p + s) - Mll\l5,p + q, k, s, k + p + q + s) 
or Q'^toV) ip,q,k,s,p + q + k + s) = (a4]}) {p, k, s, k + p + s) - (a4])) {p + q, k, s, k + p + q + s) 

V / aaaaa \ J aaar \ J raaa 
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B. Full Vertices 



We verify that the ward identities derived above hold for the fuU 3-point and 4-point vertex functions. 



1. 2-point vertex function 

We give the results for the 2-point vertex function shown in Fig © . These expressions will be needed to verify the 
ward identities for the 3-point vertex function. 

1 



1+P p 



FIG. 6: The 2-point vertex function 



E(2,p) = ^ 




f 

(2^)4 


a{l)ril+p)T{2,pJ+p) {T{4,l+p,p) 


+ r{3,l+p,p)NF{l- 




-T{2,1+p,p)Nb{1)) 


S(3,p) = - 
S(4,p) = - 




(2^)4 

f dH 

(2^)4 


a{l+p)ril)T{5J+p,p) iT{7,p,l+p) 


-T{3,p,l+p)NFil- 


I-P)- 


^T{5,pJ+p)Nb{1)) 



[a{l+ p)r{l)T{5, l+p,p) (r(8, p,l + p) - r{i,p,l + p)Nf [l + p) + [T{&,p,l + p) -T{2,p,l + p)Nf{1 +p))Nb (/)) 
+ a{l)r{l+p)V(2,p,l+p) {T{%,l+p,p) + T{7,l+p,p)NF{l+p)-{T{p,l+p,p)+V{b,l+p,p)NF{l+p))NB{l))] 



C. 3-point vertex function 

First we calculate all seven components of the three graphs shown in Fig lO that contribute to the 3-point vertex 
function. We give the result for one example: r[2,p,p-|- q]. 




p+q 




p+q 



triangle 



leaf-A 



leaf-B 



FIG. 7: The 3-point vertex function 



i f dH 

rtriangio(2,p,p + g) = - j a{l)r{l + p)T{2,p,l + p) 

{a{l + s)r(3, / + s, s)(r(6, l+p,l + s)+ r(5, l+p,l + s)Nf{1 +p)- T{2, l+p,l + s)Nf{1 + s)) 
+ r{l + s)T{2, l+p,l + s)(r(4, l + s,s) + r(3, / + s, s)Nf{1 + s) - r(2, 1 + s, s)Nb{1))) ] 



22 



ricafA(2,p,p + g) = - 
[ a{l)r{l + s)M 

rieafB(2,p,p + g) 



dH 



2 J (27r)4 

a{l)r{l + s)M{2,p, q,l,l + s){T{A, l + s,s) + r(3, 1 + s, s)Nf{1 + s) - T{2, 1 + s, s)Nb{1)) ] 

1 r dH 

2 J (2^ 

[a{l)r{l + p)V{2,p,l + p){M{<6J + p,q,~l,s) + M{h,l + q, -I, s)Np{l + p) - M{2,1 + p, q, -I, s)Nb{1)) ] 

Next, we verify the ward identities for the seven vertex functions. We give detailed resuhs for one example: 
T[2,p,p + q]. Starting from (|86|l and contracting with Q we obtain: 



where 



m 

x[2] 

y[2] 



dH 



Q-rtriangic(2,p,p + 9) = x[2]+y[2] 

Q-rieafA(2,p,p + <z) = -X[2]+li[2] 
Q-rieafB(2,p,p + g)= -y[2]+y[2] 



a{l)r{l+p + q)T{2,p,l+p) 



(86) 



2 J (27r)4 

(r(4, l+p + q,p + q)+ r(3, l+p + q,p + q)NF{l + P + q) - T{2, l+p + q,p + qjNgil)) 
= J-^[ami+p)Ti2,p,l+p) 
(r(4, l+p + q,p + q)+ r(3, l+p + q,p + q)NF{l + p) - T{2, l+p + q,p + g)A^s(O) ] 

j ^^[a{l)r{l+p + q)T{2,p + q,l+p + q) 
(r(4, l+p + q,p + q)+ r(3, l+p + q,p + q)NF{l + P + q) - T{2, l+p + q,p + q)NB{l)) 

i r dH 

= "27 (2^ 



(87) 



ail)r{l+p)Ti2,pJ+p) 



(r(4, i + p,p) + r(3, 1 + p,p)Nf{i + p)- r(2, i + p,p)Nb{i)) ] 

Comparing with (|85() we obtain 

Q ■ r(2,p,p+ g) = Q ■ (^rtrianglc(2,P,p + q) + ricafA(2, p, p + (?) + TicafB (2, p, p + q)) = S(2,p) - E(2,p + g) 

which agrees with (|83(l . The results for all components agree with the results listed in (|82|l . 

i. 4-poi'n.t vertex function 

We verify the ward identity for the 15 4-point vertex functions. We give detailed results for one example: 
M{2,p,q,k,p + q + k). There are five types of diagrams to consider. They are shown in Figs. lO, H1U|I . 

(EH), dni). 



p+q+k 
— * — 



p+q+k 



FIG. 8: The box and crossed-box diagrams 




FIG. 9: The leaftail diagrams 



p+q+k p+q+k 



/ \ P./ \ 

q' k q' k 

jellyfish-l jellyfish-2 



FIG. 10: The jellyfish diagrams 




P 1 p+q+k P 1 p+q+k 

seagull- 1 seagull-2 



FIG. 11: The seagull diagrams 




FIG. 12: The polywog diagram 
Contracting with Q and using H86|l and (|87|l we obtain: 

J- seagull— 

i(2,p, g, k, u) = a[2] -f A[2] ; rscaguii-2(2,_p, g, u) = /3[2] + _B[2] 

-'^ leaftail— 1 

(2,p,q,k,u) = 7[2] + C[2] ; 

leaftail— 3 

(2,p, g, = (5[2] +L)[2] 

^ jellyfish— 

i(2,P, <7, = e[2] + E[2] ; 

-'^ jellyfish— 2(2, p, q,k,u) = 0[2] +F[2] 

-L polywog 

{2,p,q,k,u) = G[2]+H[2] 

ricaftail-2(2,p, 9, fc, u) = -B[2] - F[2] ; rieaftail-4(2,p, q, k, u) = -A[2] - E[2] 
rbox(2,p, q, k, u) = -G[2] - C[2] ; r„.o,,ed-box(2,p, g, k, u) = -H[2] - D[2] 

where 

a[2] = -ricafB(2,p + (7,M) 
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/3[2] =rieafA(2,p,t) 

7[2] = ~rtnanglc(2,p + g,w) 

6[2] = rtrianglo(2,p, i) 

e[2] = ricafB(2,p,i) 

(j)[2] = -Tica{A{2,p + q,u) 

^ ^ 2 7 (2^)4 

[ a(/)r(/ + s)r(2,p,Z +p)(A/(6,; + s,A:, -?,it) + M{2,,1 + s, k, -1,u)Nf{1 + s) - M{1,1 + s,k,-l,u)NB{l)) ] 

^[2] = ^ / T^rir [ aUHl + t)M{2,p, k,l,l + t){T{4, l + u,u) + r(3, 1 + u, u)NFil + t) - r(2, 1 + u, u)Nb{1)) 
2 J [2Trp 



{a{l + u)T{3, 1 + u, u){T{6, l + s,l + u) + T{5, l + s,l + u)NFil + s) - r(2, l + s,l + u)Nf{1 + u)) 
+ r{l + u)r(2, 1 + 3,1 + u)(r(4, l + u,u) + r(3, 1 + u, u)Nf{1 + u)- T{2, I + u, u)Nb{1))) ] 

^[2]-^ / ^[<l)ril+p)ri2,p,l+p) 

{a{l + t)r(3, 1 + u, u)(r(6, l+p,l + t)+ r(5, l+p,l + t)NF{l +p)- r(2, l+p,l + t)NF{l + t)) 
r{l + <)r(2, l+p,l + i)(r(4, l + u,u) + r(3, 1 + u, u)Nf{1 + t)- r(2, 1 + u, u)Nb{1))) ] 
i f dH 



[ a{l)r{l + p)V{2,p, I + p) (M (6, l + s,k, u) + M (5, l + s,k, u)Nf{1 + p) - M (2, l + s,k, -I, u)Nb{1)) ] 
F[2] = -\j [ + u)M{2){p, k,l,l + t){T{4, l + u,u) + r(3, 1 + u, m)7Vf(^ + u) - r(2, 1 + u, u)NBil)) ] 



G[2] = - j ^[a{l)r{l+p)T{2,p,l+p) 



{a{l + u)T{3, 1 + u, w)(r(6, l + s,l + u) + T{5, l + s,l + u)Nf{1+p) - T{2, l + s,l + u)Nf{1 + u)) 
+ r{l + u)V(2, l + s,l + u)(r(4, l + u,u) + T{3, 1 + u, u)Nf{1 + u) - T{2, 1 + u, u)Nb{1))) ] 

i r dH 

2j (2^4 

{a{l + ?i)r(3, 1 + u, u){T{6, l+p,l + t)+ r(5, l+p,l + t)NF{l + p) - T{2, l+p,l + t)NF{l + u)) 
+ r{l + u)T{2, l+p,l + <)(r(4, l + u,u) + r(3, 1 + u, u)Nf{1 + u) - T{2, 1 + u, u)Nb{1))) ] • 



^[2]--:^ / 77^[ail)r{l+p)T{2,p,l+p) 



Using r = Ttriangio + TioafA + TioafB and combining we have 

■ M^,{2,p, q, k, u) = r,(2,p, t) - T,{2,p + q, u) 
which agrees with (|84|) . The results for all components agree with the results listed in (|82|l . 



X. CONCLUSIONS 



Calculations in real time statistical field theory are complicated by the extra indices that result from the dou- 
bling of degrees of freedom. Because of this technical problem, many people avoid the real time formalism of finite 
temperature field theory in spite of its significant advantages, as compared with the imaginary time formalism. 
Two of the major advantages of working in real time are the fact that analytic continuations are not necessary, 
and that it is easy to generalize to non-equilibrium situations. In this paper we have made a contribution to- 
wards reducing the technical difficulties associated with the real time formulation of statistical field theory. We 
have written a Mathematica program that performs the contractions over the tensor indices that appear in real time 
statistical field theory and determines the integrand corresponding to any amplitude. The program is designed so 
that it can be used by someone with no previous experience with Mathematica. It is available on the internet at 
www.brandonu.ca/physics/fugleberg/Research/Dick.html. It can be used in the 1-2, Keldysh or RA basis, and it can 
do calculations in or out of equilibrium. 
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We have used the program to calculate the QED ward identity for the 3-point function (2" — 1 = 7 components), 
the 4-point function for two fermions and two photons ((2" — 1 = 15 components), and the 5-point function for two 
fermions and three photons (2" — 1 = 31 components). Some of these identities have appeared previously in the 
literature, but the complete set of identities has not previously been published. The calculation therefore serves two 
purposes: it provides a check of the program, and it produces useful new information. We give a table that lists the 
results for the ward identities (Eqn. (|82|l ). In addition, we give a simple general expression for the KMS conditions 
between n-point functions and vertex functions, in both the Keldysh and RA bases (Eqns. H6U|I and ()56|l '). 
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